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1 Introduction 



In this paper we consider the 4D BF theory, aiming at the study of the 
dynamics induced on a planar 3D boundary. The boundary is treated ac- 
cording to the Symanzik's method [1] which strictly observes the Quantum 
Field Theory (QFT) rules of locality and power counting, together with the 
idea of separability which is a request concerning the propagators, i.e. 

A (j}1(j}2 (xi,X2) = ifxi-x 2 <0, (1.1) 

where <fti(xi) and ^2(^2) are an Y couple of fields of the theory. In a sense, the 
constraint ([Lip represents an operative definition of a QFT with boundary: 
the boundary separates the world into two parts in such a way that nothing 
propagates from one side to the other. 

This simple idea has been very fruitful for the analysis of the physics on a 
planar boundary in several circumstances. The original Symanzik's motiva- 
tion was the study of the Casimir effect [T], but later on it has been exploited 
in Topological Quantum Field Theories (TQFTs), which, as it is well known, 
do not have local dynamics but their observables are globally defined and 
deal mostly with geometrical properties of the manifolds they are built on 
[21 [3]. A boundary breaks Lorentz invariance and also the topological char- 
acter of the theory is lost; for this reason only in presence of a boundary a 
local dynamics might appear in TQFTs. This is the case for the 3D Chern- 
Simons (CS) theory [H El E] which, on the boundary, displays conserved 
chiral currents, which are the insertions of the fields on the two sides of the 
boundary obeying a Kac-Moody (KM) algebra whose central charge depends 
on the CS coupling constant. Similarly, the 3D BF theory presents, on the 
boundary, an algebraic structure carried by two sets of conserved chiral cur- 
rents [51 |9]. The BF algebraic structure is more complex than the CS case, 
but we shall come to this point later. 

However, in recent years, it is in the field of condensed matter physics that 
TQFTs with boundary has received much attention. This is mainly due 
to the discovery in the early 1980s of a new state of matter which cannot 
be described in terms of a symmetry breaking mechanism but in terms of 
topological order. In contrast to the Quantum Hall Effect (QHE), where the 
magnetic field breaks time-reversal (T) symmetry, a new class of T invariant 
systems, called topological insulators (TI), has been predicted |10| 111] and 
experimentally observed |12| in 3D, leading to the Quantum Spin Hall Effect 
(QSHE). At the boundary of these systems, one has helical states, namely 
electrons with opposite spin propagating in opposite directions [13j. The low 
energy sector of these materials is well described in terms of the CS model 
and the BF theory. Moreover, it is important to study these models in the 
presence of a boundary in order to analyze the dynamics of the edge states. 
In fact, it is well known that the edge dynamics of the QHE is successfully 
described by the abelian CS theory with a boundary both in the integer and 
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in the fractional regimes [14J. In the case of 3D TI an abelian doubled CS 
which is equivalent to a 3D BF theory with cosmological constant [8j, 
has been introduced [U [15]. Moreover, the authors of |15| argued that the 
abelian 4D BF theory with a boundary could describe some features of TI 
in 3D, in particular their edge fermionic degrees of freedom. 

The paper is organized as follows: in Section 2 the bulk 4D BF theory is 
described. The field equations and the Ward identities describing the resid- 
ual gauge invariance typical of the axial gauge are derived. In Section 3 
the boundary is introduced following the Symanzik's approach. The most 
general boundary action and the corresponding boundary conditions on the 
fields are written, together with the Ward identities of the theory, modified 
by the boundary. In Section 4 the boundary algebra is computed, which gives 
rise to a surprising electromagnetic structure on the boundary, which is to 
be understood both from the geometrical (x% = 0) and the field theoretical 
(mass-shell) point of view. An "electric" scalar and a "magnetic" vector po- 
tential are identified, by means of which we write a 3D boundary Lagrangian. 
A remarkable "duality" relation between the electromagnetic potentials is ob- 
tained. All the results are summarized and discussed in detail in Section 5, 
and the propagators of the full theory, including the boundary, are explicitly 
given in the Appendix. 

2 The classical theory 

In the abelian case the action of the four-dimensional BF model O [16] , which 
describes the interaction between the two- form B^ v and the gauge field A^, 
is given by: 



where = d^A u — d u A^. It is well known that the BF theories, in any 
spacetime dimensions, both in the abelian and non-abelian case, do not 
depend on any coupling constant. Here, k is a constant which we have 
introduced in order to distinguish the boundary terms from the bulk terms; 
it can be eventually put equal to one at the end of the computation. 
The action (|2.ip is invariant under the symmetries: 




(2.1) 



6® A, 
6^B, 



jus — 








(2.2) 



and 



5 {2) A a = 




(2.3) 
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where 6 and are local parameters. 

We remind that the action (I2.ip is not the most general one compatible with 
the symmetries (j2.2|) and (|2.3p . Indeed, a Maxwell term J c! iF^F" could 
be added, coupled to an additional parameter. In this paper, we consider 
the action (|2.ip alone, because we are interested in the 3D dynamics on the 
edge of a TQFT. We are allowed to do that, if we think of the action (|2.ip 
as the abelian limit of the non-abelian 4D BF theory |TT], which, as any 
other TQFT, is protected from the occurrence of non-topological terms by 
an additional symmetry, called "vector supersymmetry" in |18j . Nevertheless, 
the non-abelian case, which is much richer from the fields theoretical point of 
view |19[ [20] , and the addition of a Maxwell term are interesting extensions 
which we shall present elsewhere [21J. 

Lorentz invariance will be broken by the introduction of a planar boundary. 
Consequently, a convenient choice for the gauge conditions on the two fields 
involved, is the axial one: 



where latin letters run over 0,1,2. The axial choice is implemented by adding 
to the action (|2.ip the gauge fixing term 



where b and d 1 are respectively the Lagrange multipliers for the fields ^3 and 
B{ 3 . As usual, in the abelian case the ghost fields are decoupled from the 
other fields. One of the main differences with the non-abelian case, is the 
structure of the gauge fixing term, which, because of the reducible symmetry 
(12. 3p . involves ghosts for ghosts, and therefore is highly non trivial |19[ 120] . 
Summarizing, the classical action is given by 



r c [J^] = / dtxiKe^pdiAjBks + (diA 3 - d 3 Ai)B jk ] + bA 3 + ^B a + 



where J$ are the external sources coupled to the quantum fields <3?. From 
(|2.6p . we get the bulk field equations: 




(2.5) 



+ Jg.Bij + 2J§ a B i3 + J\.Ai + J A3 A 3 + J b b + J l di di} 



(2.6) 



J\.+Ke i ^[2d j B k3 + d 3 B jk ] = 
J% j +Ke i i k (d k A 3 -d 3 A k ) = 
J 3 + b - Ke ijk diB jk = 
2J|. 3 +d i + 2 K e ijk d j A k = 
A 3 + J b = 



(2.7) 



B i3 + 4 = 0. 
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It is well known |22| . that the axial gauge is not a complete gauge fixing. A 
residual gauge invariance remains on the plane X3 = 0, which is described 
by two (one for each gauge symmetry (I2.2p and (|2.3[) ) local Ward identities: 

w(x)r c [j^} = diJ% + d 3 j% + d 3 —^ = o, (2.8) 

W\x)Y c m = + dsJ^s + ia 3 ^f = 0. (2.9) 

In the following Table 1, we list the canonical mass dimensions of the fields 
of the theory: 







Bfj. u 


b 


d l 


Dim 


1 


2 


3 


2 



Table 1: Canonical mass dimensions of the quantum fields. 



3 The boundary 

To introduce a planar boundary in the theory, we adopt the Symanzik's 
method [I], which basically consists in writing the most general boundary 
Lagrangian, according to the general QFT principles of locality and power 
counting, and then computing, for the modified theory, the propagators, on 
which the constraint of "separability" is imposed. This corresponds to asking 
that the propagators between points on opposite sides of the boundary x = 0, 
vanish: 

A^ 1(f)2 (xi,X2) = if xi ■ x 2 < 0, (3.1) 

where <j)i(xi) and (^2(^2) are two generic fields of the theory. No other as- 
sumption is required. In particular, as we shall see, no boundary condition 
is imposed. Rather, we shall find out which are the most general ones com- 
patible with the presence of a separating boundary 

The most general boundary Lagrangian must respect locality power count- 
ing, and covariance in the plane X3 = 0. An additional, more subtle, con- 
straint comes from the fact that the field equations (|2,7p appear to be first 
order differential equations. This originates from the BF action (|2.ip . which 
contains only first derivatives. Hence, in order to be solved, only one bound- 
ary condition on any field is needed, and no more than one is allowed. This 
obvious consideration translates into the request that the boundary term in 
the action should not contain differentiated fields (with respect to X3). In 
the present case, only a term proportional to 5{x 3 )A l d 3 Ai(x) would matter, 
resulting in a boundary condition on the differentiated fields &3Ai(x)\ _ Q , 
in addition to the one on the undifferentiated fields Ai(x)\ X3=0 , which is al- 
ready present. This would be incompatible with the first order differential 
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equations (|2.7p . Hence, such a term will not be considered. Two more re- 
marks on this point: the first on the fact that this is peculiar of the 4D (or 
higher) since in the lower dimensional cases (CS and 3D BF) such a term 
violates the power counting constraint. The second concerns the possible 
presence in the bulk of a Maxwell term, which, being of second order in the 
derivatives, would spoil the previous argument, and the only possible term 
we are talking of, namely 5(x 3 )A l d 3 Ai(x), should be in that case included. 
Summarizing, the most general boundary Lagrangian is: 

Cbd = S(x 3 ) axAiB 1 + a 2 jA i A i + a 3 b + ^e iik diAjA k + a^A 1 , (3.2) 

where B l = ^ k Bj k (and correspondingly J^ ih = e l ^ k J^ k ), and a a ,a = 
1, . . . , 5 are constant parameters to be determined. In order to have all the 
a-parameters massless, a 2 has been given an explicit m-mass dependence, 
which will turn out to be important in the following. 

The separability condition (|3.ip allows for a considerable simplification, by 
considering only one of the two sides of the boundary. In fact, basically, the 
constraint (I3.ip means that the opposite sides of the boundary are completely 
decoupled. Hence, we can consider just one side of the boundary, say the '+' 
side, forgetting about the opposite side, which can be obtained from the '+' 
side by parity. 

Having this in mind, the boundary Lagrangian Cm (|3.2p modifies the bulk 
field equations (|2.7p as follows: 

J\. + 2Ke ijk djB k3 + ad 3 B i = -5{x 3 )[ ai B i+ + a 2 mA i+ + a 3 b + 
+ a^ k {d j A k )+ + a 5 d i+ ] 

& k J §k + ^ k {d k A 3 - d 3 A k ) = - ai 5(x 3 yi k A+ 

J 3 + b - Kdi& = (3.3) 

2jg 3 + cf + 2Ke ijk djA k = 

A 3 + J b = -5(x 3 )a 3 

B i3 + 4 = -5(x 3 )a 5 A l+ , 

where the apex + denotes the insertions of the fields of the theory on the 
'+'-side of the boundary x 3 = 0. 

Consequently, the boundary term Cm (|3.2p breaks the local Ward identities: 
diAi + ds4 3 + ^ = -5{x 3 )[ ai diB i+ + a 2 mdiA i+ + a 5 8id i+ l (3.4) 

e^dj J §k + d 3 jg 3 + ^d 3 d l = -5(x 3 )a ie ^ fc 9^+ (3.5) 

We stress that the breaking terms at the r.h.s. of (|3.4p and (|3.5p are linear in 
the quantum fields, and hence a non-renormalization theorem ensures that 
they are present at the classical level only, and do not acquire quantum 
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corrections |23| . For this reason, linearly broken symmetries are perfectly 
allowed by general principles of QFT. 

Postulating that b{x 3 = ±00) = d l (x 3 = ±00) = 0, the Ward identities (|3.4p 
and (|3.5p can be integrated as follows: 



3.1 The boundary conditions 

When dealing with a boundary, the question arises of the boundary con- 
ditions on the fields and their derivatives. Symanzik's approach allows to 
find out the most general ones, without imposing any additional constraint. 
In fact, the separability condition ([3.ip and the boundary Lagrangian (|3.2p 
can be used to calculate the propagators everywhere, in the bulk and on its 
boundary. Studying the limit of the propagators on the different sides of the 
boundary, one can infer which are the possible boundary conditions on the 
fields. 

This program has been carried out for the 4D scalar QFT [TJ [3], for the 
CS theory [5j [6] and for the 3D BF model [8]. But it might happen that 
the explicit calculation of the propagators of the whole theory (bulk and 
boundary) may become very difficult (this is the case for the 3D Maxwell- 
Chern-Simons theory with boundary |24j . for instance). Nonetheless, even in 
that case it is possible to find out the boundary conditions and the physics on 
the boundary (typically the existence of conserved currents and the algebra 
they form), avoiding the explicit calculation of the propagators. 
In this paper, we are mainly interested in the dynamics of the 3D boundary of 
the 4D BF theory, so we proceed directly to determine which are the possible 
boundary conditions which characterize the physics on the boundary. The 
complete propagators of the full theory are given in the Appendix. 
In order to find out the most general boundary conditions, we integrate the 
broken equations of motion (|3.3p in an infinitesimal interval around X3 = 
|24j . Because of the separability condition (|3.ip . we get the following 
algebraic system, involving quantities lying on the '+' side of the boundary 
(the opposite one being obtained by parity), i.e. the fields on that side of 
the boundary and the a-parameters: 




(3.6) 




(3.7) 



{K + a x )B i+ = - 
{K-a 1 )Af = 
a 3 = 
a 5 A i+ = 0. 



CLlVflA 



i+ 



a 3 b + - a A e ijk djA+ - a 5 d' 



(3.8) 
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The solutions of the system fj3.8 j) are acceptable if the r.h.s. of the Ward 
identities fj3-6|) and (|3.7p do not vanish. Indeed, in such a case, differentiating 

(or ()3.7p ) with respect to J a (or J^). we obtain the 



the Ward identity f[376|) 
inconsistency 







(3.9) 



The interesting fact we learn from this simple remark, is that the boundary 
Lagrangian (|3.2p acts as a kind of gauge-fixing for the residual gauge 
invariance, in the sense that its presence is necessary (and sufficient) to 
calculate the propagators on the boundary. 

We see that, if we ask that the boundary term of the Ward identity (|3.7p 
does not vanish, we must impose the condition A l+ 7^ 0. Consequently, it 
must be: 



a 5 = 
ai = k. 



(3.10) 



Having done that, the system (|3.8p reduces to the single equation: 

2kB 1+ = -a 2 mA i+ - a 4 e ijk djA+, (3.11) 
which has four different algebraic solutions, listed in Table 2. 





a2 


04 


A i+ 




1 


*o 







«2m 41+ 
~2kT 


2 













3 









_04 ijkf). A + 


4 


^0 


^0 




-£(a 2 mA*+ + a 4 e« fc a 7 -A+) 



Table 2: Algebraic solutions of the equation (|3.1ip 



Now we notice that solutions 2 and 3 lead to ill-defined Ward identities (the 
r.h.s. of (|3.6p vanishes), and, for this reason, they are not acceptable. 

We are then left with two sets of boundary conditions which satisfy the 
system (|3.1ip : 





(12 


04 


A i+ 


B i+ 


I 







7^0 


aim Ai+ 


II 


7^0 


7^0 


7^0 


-±(a 2 mA l+ + a^ k djA+) 



Table 3: Acceptable solutions of the equation (|3.1ip 



In both cases the parameter a 2 must be different from zero: the presence 
of the massive term m is necessary in order to make the theory consistent. 
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In other words, it is necessary that the boundary Lagrangian is not scale- 
invariant. This fact is very important for the study of the physics on the 
boundary, as we shall see later. 



4 The algebra and the physics on the boundary 

In this section we derive the algebra of local observables which is generated 
on the boundary, due to the residual gauge invariance of the theory, func- 
tionally described by the Ward identities (|2.8p and ()2.9[) . 
Next, we shall argue that it is possible to describe the physics on the bound- 
ary in terms of two fields: a gauge field 0(^0 an d a scalar massless field 
A(X) (we recall that X = {x$,x\,X2), while x = (X, X3)). We shall identify 
the Lagrangian which describes the physics on the boundary by interpreting 
the boundary algebra as a set of canonical commutation relations for the 
fields Q % and A. 

4.1 The boundary algebra 

Despite the fact that the solutions listed in Table 3 appear to be different, 
and depending on free parameters, the broken Ward identities for both of 



We therefore remark that, as a matter of fact, a unique solution exists, which 
does not depends on any free parameter. We recall that the constant k was 
introduced in order to keep trace of the bulk dependence, but it is not a 
true coupling constant. We can therefore freely put n = 1 in what follows. 
Evaluating the previous relations at vanishing sources, i.e. on the mass shell, 
we find that: 



We now differentiate fj4. 1|) with respect to J l A (x'), with x' lying on the '+' 
side of the boundary X3 = 0, obtaining: 



them are: 




(4.1) 



(4.2) 



diB i+ = 
e ljk djA+ = 0. 



(4.3) 
(4.4) 



W^iX' -X) = di(^ Ai&i {x',x)) 



.0+ 



(4.5) 



Next, we express the propagator in (|4.5j) as follows: 



(A AiSi (x',x))^ =x/=o+ = e(t-t')(B i+ (X)A+(X'))+0(t> 



t)(A+(X')B*+(X)). 

(4.6) 
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Substituting the previous identity in (|4.5p . we find: 

Sjd^iX' -X) = S(t - t')([B°+(X),A l (X')])+ 

+ (9(t - t')(d i B i +(X)A+(X')) + 0(t' - t){Ai(X')diB l+ (X))). 

(4.7) 

Remembering (|4,3p . we obtain: 

5(t-t')[B 0+ (X),A l (X')} = 5jd i 5^(X' -X). (4.8) 

For Z = 1,2, it is possible to factorize 5(t — t'), finding: 

[B°+(X),A a (X%= t , =d a 5^{X' - X), (4.9) 

where a denotes the indices 1,2. From now on, greek letters will denote 
spatial coordinates on the plane X3 = 0. 

Next, differentiating (|4.ip with respect to J^(x'), with x' lying on the '+' 
side of the boundary, we get: 

a t (A^(x',x)) ; =0. (4.10) 

Following the same reasoning which led to (|4.9p . we obtain: 

[B 0+ (X),B+(X% =t ,=0. (4.11) 

In particular, if I = the previous commutation relation become: 

[B 0+ (X),B+(X% =t , = 0. (4.12) 



Let us now consider (|4.2p . The differentiation of this identity with respect 
to J A i(x'), (with x' 3 = + ), leads to: 

e^(A AlAk (x',x)) = 0. (4.13) 

Taking into account f)4.4[) . the previous identity yields the following commu- 
tation relation: 

[At(X),A^(X% =t ,=0. (4.14) 
Next, differentiating (|4.2p with respect to J~(a/), we find: 

(d j d l k -d k 5 l j )5^(X'-X) = 

-^(A^ fc (*',x)) +d k (A &A (x',xj) ; 

V / X3=x 3 =0+ V 3 / X3=Xg=0+ 

(4.15) 
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which does not provide new commutation relations. 

In conclusion, the commutation relations (|4.9p . (|4.12p and (|4.14|) form the 
following algebra of local boundary observables: 

[B 0+ (X),A a (X')] t=t , = d a 5^(X' - X) 

[B 0+ (X),B+(X% =t ,=0 (4.16) 

[A+(X),AP+(X%= t ,=Q, 

which will be discussed in detail in the last section of this paper, together 
with the other results. 



4.2 The physics on the boundary 

Let us consider again the broken Ward identities (|4.1|> and (|4.2p . As we 
said, they describe the residual gauge invariance of the 4D BF theory on the 
planar boundary X3 = 0. Going on the mass shell, i.e. at vanishing external 
sources, we find the equations (|4.3p and (|4.4p for the fields on the boundary, 
for which the conditions listed in Table 3 hold. The equations (|4.3p and (|4.4p 
are easily recognized as the electromagnetic Maxwell equations for an electric 
(-R- B l+ ) and magnetic (0 ^4^) field. In other words, the 4D BF theory, 
which, as any other TQFT does not have local observables and has vanishing 
Hamiltonian, when dimensionally reduced on a 3D planar boundary, acquires 
a rich, physical, electromagnetic structure. 

But we can push this further. The conditions (j4.3[) and (|4.4p allow us to 
express the fields B l+ and A l+ in terms of the potentials A and Q: 

d r B i+ = => B i+ = e^djCk 
^%A\ = A+ = d k A, 

where A(A) and C*(A) have canonical dimensions zero and one, respectively. 
The fields Af and B l+ are left invariant by translational and gauge trans- 
formations of the potentials as follows: 

5 A = c 

(4.18) 

where c is a constant and 0(X) is a local parameter. 

Let us now consider the boundary condition / in Table 3. With a suitable 
choice of the free parameter 02 , we can rewrite this condition in terms of the 
fields Q l and A: 

e ijk djCk = md l K. (4.19) 

The massive parameter m in equation (|4.19p allows a rescaling of the fields 
Q and A as follows: 

A 

A —;= 

m (4.20) 
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So, thanks to the massive parameter m, the rescaled fields can be given 
the standard canonical dimensions of a gauge field and of a scalar field in 
three space-time dimensions ([£*] = [A] = ^). Therefore, the equation (I4.19P 
becomes: 

e ij %Ck = d l A, (4.21) 

which is exactly the duality relation between a scalar field and a gauge field 
which is required to construct massless fermionic fields in three dimensions 
via the tomographic representation [25]. This could be interpreted as the 
sign that the actual degrees of freedom of the 3D theory obtained on the 
boundary are fermionic rather than bosonic. We shall come back to this 
point in the conclusive Section 5. We now consider the solution II in Table 
3: 

B i+ = -(a 2 mA i+ + a 4 e ijk d j A+). (4.22) 

It is evident that, if we evaluate the previous condition on the mass-shell, the 
term proportional to vanishes due to the condition (|4.4|) . and the previous 
equation is equivalent to the boundary condition /: 

B i+ = -a 2 mA i+ . (4.23) 

In other words, the duality condition (I4.2ip always holds, and the 3D physics 
we are discussing here is therefore uniquely determined. 
We are now able to find a 3D Lagrangian for the fields and A which 
describes the physics on the boundary and which is compatible with the du- 
ality condition (|4.2ip . In what follows, we shall interpret the algebra (|4.16p 
as a set of canonical commutation relations for the fields Q and A, and we 
shall find the corresponding Lagrangian, doing the contrary of what is com- 
monly done, which is to find the canonical variables and their commutation 
relations from a given Lagrangian. 
Let us consider the equation (j4.8|) with / = 0: 

5(t - t')[B 0+ (X),A (X')] =5'(t- t')5^(X - X'). (4.24) 

Writing this identity in terms of the fields Q and A, we obtain: 

5(t - t')d' [e^d a ^(X),A(X')} = 5'(t - t')8^(X - X'), (4.25) 

where we have factorized the operator d' Q on the right hand side since it acts 
only on the field A. It is easy to see that 5(t — t')d' Q = —5'(t — t') and, con- 
sequently, we can factor out the 5'(t — t'), finding the following commutation 
relation: 

[A(X'),e a Pd a (p(X)] t=t , =5^{X-X'). (4.26) 
Consider then the first commutation relation in (|4.16p : 

[B 0+ (X),A a (X')] t=t , = d a 8^(X'-X). (4.27) 
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If we express the previous identity in terms of the fields A and we find: 



d a [e^Cl3(X),d' y A(X% =t , = 5^d v 5^(X' - X), (4.28) 

which yields: 

[e a %(X),<9;A(X')W = 5^ 2 \X'-X). (4.29) 

We are now ready to construct the Lagrangian. The commutation relations 
(|4~26|) and (pi~2"9|) allow us to interpret the fields II (A) = e a/3 <9 a C/3 and U^ a = 
d a A as the conjugate momenta of the fields A and £° = e a ^Cp respectively. 
With these assumptions, the Lagrangian of the system C = Yl n<& — H , 
where H is the Hamiltonian of the system, is given by: 

C = e^d a (pd t A + daA^dtCp - (e a Pd a (p) 2 - (d a A) 2 , (4.30) 

which is equivalent to the Lagrangian postulated in [15J for the study of 
the topological insulators. We stress that the Lagrangian (|4.30p is the most 
general one compatible with power counting and respecting the symmetries 
Furthermore, the coefficients of the terms appearing in (|4,30p are 



fixed by making the field equations of motion compatible with the duality 
relation (|4.2ip . Moreover, if we omit thekinetic term, C is equivalent to the 
Lagrangianconsidered in |26] to study the edge states of the 4D BF theory. 



5 Summary and discussion 

The main results presented in this paper are 
1) boundary as gauge fixing 

According to the Symanzik's approach, the boundary conditions on the fields 
are not imposed, but, rather, are derived from the form of the propagators. 
Now, this is not always feasible, and almost always quite difficult (we are 
talking about computing the propagators of the theory, including the bound- 
ary and satisfying the separability condition (jl.ip ). In a previous paper |24| 
we already succeeded in finding out the boundary conditions (which are the 
starting point for the study of the physics on the boundary) without com- 
puting explicitly the propagators. At the end, one is faced with a nonlinear 
algebraic system whose unknowns are the parameters on which the bound- 
ary Lagrangian depends, and the fields (and its derivatives) on the boundary. 
Most of these solutions are inconsistent, or unacceptable for some reasons, 
and in |24] these unphysical solutions were ruled out one by one. Here, we 
found a nice, general criterion to get the same result: the r.h.s. of the Ward 

1 We thank the referee for suggesting us to clarify the issue of the uniqueness of (|4.3U|) 
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identities (|5.ip and ()5.2p . i.e. the linear breaking due to the boundary term 
in the action, must always be different from zero. Otherwise, the propaga- 
tors cannot be defined. The nice interpretation of this statement, is that the 
boundary term in the action plays the role of a gauge fixing of the residual 
gauge invariance on the boundary. This observation leads immediately to 
the solutions / and listed in Table 3. 

2) Ward identities in presence of a boundary 

/+oo 
dx 3 diJ% = 8iB i+ (5.1) 
-oo 

/+oo 
dx 3 e i ^d j J §k = -e l3k djA+. (5.2) 
-oo 

Quite remarkably, the apparently distinct solutions / and II of Table 3 phys- 
ically coincide, since they lead to the same Ward identities on the boundary. 
This is the first evidence of the striking electromagnetic structure which 
determines the physics on the boundary, as we shall discuss shortly. In ad- 
dition, despite the fact that the solutions depend on free parameters, when 
put into the Ward identities (|5.ip and (|5.2p . which contain all the physical 
information, these disappear. The separability condition isolates a unique 
dynamics on the boundary, without any dependence on free parameters. 



3) electromagnetism on the boundary 



diB 1 



=> B t+ = e ijk djCk (5.3) 
=> A+ = d k A, (5.4) 

On the boundary x 3 = 0, and on the mass shell = = (we stress 

this double constraint defining the boundary), the 4D topological BF the- 
ory displays Maxwell equations for an electric field and a magnetic field, 
to be identified with the boundary insertions B t+ and Af , respectively. 
This is a direct consequence of the result 2). Consequently, two potentials 
can be introduced: an electric scalar potential A(X) and a magnetic vec- 
tor potential £*(X), depending on the 3D coordinates on the plane x% = 0: 
X = (x ,x 1 ,x 2 )- 



4) duality 

e^djCk = d'A. (5.5) 

The solutions of Table 3, i.e. the possible boundary conditions on the fields, 
translates in the "duality" condition between the potentials (|5.5p . This con- 
firms the fact that the dynamics on the boundary is uniquely determined by 
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the Ward identities (|5.ip and <|5.2[) . We find here, in a well defined field the- 
oretical framework, a strong motivation for a relation which is known since a 
long time [25j, where this duality (or "tomographic") relation was introduced 
to give a Bose description of fermions in 3D. Here, this condition appears as 
the unique boundary condition on the fields A l+ = B l+ , written in terms of 
electromagnetic potentials defined by the boundary Maxwell equations (|5.3p 
and (|5.4|) . This strongly suggest that the actual degrees of freedom of the 
dimensionally reduced 3D theory are fermionic, confirming recent develop- 
ments concerning the edge states of topological insulators, which seem to be 
described in terms of fermion fields |27) . 

5) 3D boundary algebra 

[B°+(X),A a (X% =tl = d a 6®(X' - X) 
[B°+(X),B+(X% =t ,=0 (5.6) 
[A+(X),AP+(X% =t ,=0, 

On the boundary, the above algebra is found. It is formed by a vectorial, 
conserved current, whose 3D components are the insertions of the fields on 
(one side of) the boundary (B 0+ (X) and A^(X) a = 1,2, related by the 
duality-boundary condition (|4.23p ). We stress that the conservation of the 
current is obtained on the mass-shell, i.e. at vanishing external sources J^. 
This is in perfect analogy with what happens in the 3D CS and the 3D 
BF theory. In all cases, the conservation of the currents comes from the 
Ward identities of the residual gauge invariance broken by the most general 
boundary term respecting Symanzik's separability condition (|l.ip . going on 
the mass-shell, and exploiting the boundary condition previously found on 
the quantum fields. The physical interpretation of the current conserva- 
tion is different, since in the 3D CS and BF cases, it leads, thanks to the 
boundary conditions, to the chirality of the currents. In the 4D BF case 
the current conservation (|5.3p (again, together with the duality-boundary 
condition (|5.5p ). is tightly related to the electromagnetic structure and the 
consequent determination of the electromagnetic potentials. One more com- 
ment on the algebra (|5,6p is on order. The 3D BF theory shows two types 
of algebraic structures on its 2D planar boundary, whether the cosmological 
constant A, whose presence is peculiar of the 3D case, is vanishing or not 
[8]. For A / 0, one finds the direct sum of two KM algebras, and this is 
not surprising because the 3D BF theory with cosmological constant can be 
written in terms of two CS theories with opposite coupling constants, each 
of which shows a boundary KM algebra. On the contrary, the algebra found 
for vanishing cosmological constant is truly BF-like. Indeed in that case the 
BF theory, in any dimension, cannot be rephrased in terms of CS actions. A 
remarkable check of our results, is that we find exactly the same algebraic 
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structure as the one found for the 3D BF case for A = 0, whose interesting 
features and relationships with other bulk theories are discussed elsewhere 
|21j . In a different framework and language the same algebra has been found 
in |26| . written in terms of the same dynamical variables we treated in this 
paper, i.e. the edge states of 4D BF theory. 

6) canonical commutation relations and dimensional reduction 

[A(X),n {A) (X')] t=t , = 6®(X-X') (5.7) 
~_C(X),Il {0 p(X')] t=ti = 5 a ^ 2 \X-X% (5.8) 

where 11(a) = ^daCfi an d ^-(^)a = d a A are the conjugate momenta of the 
fields A and Q a = e a ^Qp respectively. The point to stress here, is that, writ- 
ten in terms of the electromagnetic potentials (I5.ip and (|5.2p , the boundary 
algebra (|5.6[) can be interpreted as a set of canonical commutation relations, 
for the canonically conjugate variables. Once realized this, it is almost im- 
mediate to write down the corresponding 3D Lagrangian, which is uniquely 
determined by our procedure. Indeed this analysis can be viewed as a sys- 
tematic way to find (D— 1) -dimensional Lagrangians out of D-dimensional 
bulk theories. It is a surprising and welcome result, that this new way of 
dimensionally reducing D-dimensional theories originates from the algebraic 
structure found on the boundary, interpreted as a set of canonical commu- 
tation relations, and which comes from the Ward identities describing the 
residual gauge invariance on the boundary and broken (by the boundary it- 
self) in the most general (and unique) way compatible with the Symanzik's 
simple criterion of separability. 



7) 3D Lagrangian 



c = ^2~n<l>-H 

= e a Pd a (pd t A + Q Ae°^C/J " {^d a Cpf - {d a kf (5.9) 

This is the 3D Lagrangian obtained on the mass-shell boundary of the 4D 
topological BF theory. It is the unique solution compatible with the QFT 
request of locality, power counting and with the Symanzik's criterion of sepa- 
rability (jl.ip . It is left invariant by gauge and translational transformations. 
It is non-covariant, and its dynamical variables (scalar and vector potentials) 
are coupled in a non-trivial way. Quite remarkably, this action, uniquely de- 
rived here by very general QFT principles, coincides with the one studied 
in |26| for the edge states of the 4D BF theory, where, the same algebraic 
origin is stressed. In a completely different theoretical framework, the action 
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(|5.9p is employed to study the surface of 4D (3+1) topological insulators |15j . 
The duality relation (|5.5p is there exploited to extract the desired fermionic 
degrees of freedom. 

There are several interesting developments of the results presented in this 
paper. The most obvious is the non-abelian extension, with particular atten- 
tion to the duality relation (|5.5|) . to the boundary algebra (|5,6p and to the 
boundary Lagrangian (|5.9|) . Moreover, it is of interest to apply our method 
to dimensionally reduce 5D bulk theories, in order to find out the resulting 
4D actions. 



A The propagators 



In this appendix we shall derive the propagators of the 4D BF model (|2.ip 
with boundary (|3.2p . taking into account the boundary conditions in Table 3. 

The separability condition (|3.ip is satisfied by propagators with the following 
form: 



(x, x') = 6(x 3 )e(x' 3 )A+ irt>2 (x, x') + 8(-x 3 )6(-x 3 ) A"^ (x, x'), (A.l) 



where A , , are the propagators on the 'db' side of the boundary. They 
are solutions of the system of equations obtained by differentiating the bulk 
equations of motion (|2.7p with respect to the sources of the fields. They 
must be compatible with the boundary equations of motion (|3.3p and with 
the Ward identities (|3.4p and (|3.5p . Since At, and A^ 1( ^ 2 are transformed 
into each other by a parity transformation, in this appendix we derive a 
solution for At^, where x 3 ,x'^ > 0. In what follows we shall omit the 
index +. 

If we differentiate the equations of motion (|2,7p with respect to the sources 
of the fields and we evaluate the expressions obtained at vanishing sources, 
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(A.2) 



we get a system of equations for the propagators of the theory: 

A A ^{x,x') = V i>{x')^b{x') 

A A3b (x,x') = -5^(x-x') 
A B i3^(x, x') = V tfj(x') ^ d l {x') 

A B , Sd[ (x,x') = -5 i l 5^\x-x') 
d 3 A Ai§t (x',x) = -5\5^(x-x') 

d ^&Bi( x '^ x ) = 
d 3 A b Bi(x',x) = 

d 3 A di§l (x',x) = 2e l \d J 5^(x-x') 

d3A AlAl (x',x) = 

d 3 A §iAi (x',x) = 5 i l S^(x-x') 

d 3 A bA ,(x',x) = -d t 8^{x' -x) 
dsA dlAl (x',x) = 
A Alb (x',x) = diA Ai 2i(x',x) 

^Bib( X ' a: ) = ^i^ L B i B i ^ X ' • r ) 

A 66 (x',x) = diA b gi(x',x) 
A dlb (x',x) = diA d[ ^(x',x) 
Aa ; * = -26^A^(x',x) 
A^ di (x',x) = -2e^A^ Afc (x',x) 
A d;di (x',x) = -2e ij ' fc a ; -A <Mfc (x , ,x) 
A Ml (x / ,x) = -2e^a i A Mfc (x',x). 

Notice that it follows directly from the gauge conditions, i.e. from the last 
two equations in (|2.7p . that the Green functions containing ^3 and Bi 3 vanish 
except A A3 b(x,x') = —S^{x — x') and A B a dl {x,x') = —5\5^(x — x') and, 
for this reason, we do not list these propagators in the following. 
The most general solution of the previous system is: 

A^ 2 (x',x) = 

( 5/(X,X') SjT cl (x,x') -2eV i ^ fc (X',X) -^T ci (x, x') \ 

-SjT C2 (x',x) ni(X,X') 2e ij l d j T C2 (x',x) A 7 ) 

-2e%%H^(X',X) 2e ij l d j T C4 (x,x > ) 4e^ ejfcpg ^Hf (X, X') 

^ ^Tc^x'^x) d l ttt(X\X) d l d i n i l (X,X')J 

(A.3) 
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The labels (p\ and 4>2 run over the set of fields {A 1 , B l , d l ,b}. T Ci (x, x') is the 
tempered distribution (e(x 3 -x' 3 ) + c i )S 3 (X -X'), E/(X, JT) and ft/(X,X') 
are generic functions of the transverse coordinates X = (xo, xi, X2), and 
Ci,i = 1, ...,4 are constant parameters. 

Let us now consider the boundary conditions / and II in Table 3. Remark- 
ably, for both these solutions the Ward identities (|3.6p and (|3.7|) take the 
following form: 

00 

dxsdiJ^i = d t B l+ 
« (A.4) 

-00 

Differentiating the equations (|A.4j) with respect to the sources J4Z (a/), J^i (x r ), J d i (x') 
and Jb(x'), we obtain eight differential equations for the propagators: 

W^iX'-X) = dJA A ^(x',x)) (A.5) 

\ 1 / £3=0 

diU Bfii {x',x)) = (A.6) 

\ 1 / Z3=0 

djA dl§i (x',x)) = (A.7) 

\ / £3=0 



^(a^(x',x)) = (A.8) 

e^%(A AiAk (x',x)) = (A.9) 

e^ k djS kl S^(X' -X) = -e^ k dJA, A (x',x)) (A.10) 

\ 1 K / £3=0 

e^(A, ;Afc (x',x))^ =o = (A.ll) 
j1 ] >d j (A bAh ( X f,xj) = 0. (A.12) 

Substituting the propagators (|A.3P in the above system of differential equa- 
tions, we get the following constraints on the parameters c\ and c 2 : 

(A.13) 

c 2 = 0. 

Additional constraints come from the request that the Lagrange multipliers 
b(x) and d l {x) vanish at £3 — > 00. Consequently, the propagators involving 
these fields must satisfy: 

lim 2e ij ,djT C4 (x,x') = 
*3->+oc (A _ M) 

lim diT c Jx ,x) = 0, 

X3— > + 00 
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which yield: 

c 3 = -1 

(A.15) 

C4 = 0. 

Moreover, the propagators A AiA i(x' ,x) and §i(x' ,x) must be symmetric 
for the exchange {x,i} <-> {x',l}. As a consequence, taking into account 
(TO)) and (|A~TT|) . we obtain that the functions ~/(X,X') and 
take the following form: 

Ei{x,x') = d i d l 7 1 (x -x') 

J (A.16) 
nf(X,X') = eV%ef s d r <j> ks (X - X f ), 

where rj(X — X') and 4>ks(X — X') are generic functions of the transverse 
coordinates X — X' , and have canonical mass dimensions zero and two, 
respectively. In terms of these functions, the matrix of propagators (|A.3P 
finally reads: 

A AB (x',x) = 

(d l dir]{X - X') -8\T- X {x,rf) -d;T_i(x, x')\ 

-SjT (x',x) & k dje{»dr<l> ks {X-X') 2e ij l d j T (x',x) 

2e ij l d j T (x,x') 

^ d l T^i(x',x) J 

(A.17) 
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